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It has been shown how on-sheh forward scattering ampHtudes (the "Bar- 
ton expansion") and quantum mechanical path integral (QMPI) can both 



be used to compute temperature dependent effects in thermal field theory. 



^ ■ We demonstrate the equivalence of these two approaches and then apply the 



QMPI to compute the high temperature expansion for the four- point func- 
tion in QED, obtaining results consistent with those previously obtained from 

the Barton expansion. 
ll.lO.Wx 



I. INTRODUCTION 

Temperature dependent effects can be computed in thermal field theory in a variety 
of ways. In the real-time formalism, directly resorting to the Feynman diagram approach 
results in awkward sums over so-called Matsubara frequencies It is possible to perform 
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an analytic continuation in the time variable in this formalism to imaginary time, and to 
then convert the sum over Matsubara frequencies into a contour integral. The temperature 
dependent effects in one-loop processes are thereby related to an on-shell forward scattering 
amplitude weighted by a statistical factor. This is the so-called "Barton expansion" 0,0. 
An alternate approach involves regulating the one-loop generating functional using the zeta 
function [^,0 and then computing matrix elements of the form < x\e~^^\y > that arise in 
this procedure by a quantum mechanical path integral (QMPI). The QMPI is over a space 
compactfied in the temporal direction; in place of a sum over Matsubara frequencies, a sum 
over winding numbers is encountered 0. 

These two approaches to computing radiative effects in thermal field theory (the Barton 
expansion and the QMPI) are apparently quite disparate and there is no obvious connection 
between the two techniques. In this paper, we demonstrate that in fact they are equivalent. 
A calculation of the one-loop amplitude with four external photons has been computed in 
thermal QED using the Barton expansion 0; in this paper we reconsider this computation 
using the QMPI and demonstrate that the results of these approaches are compatible. 

II. THE BARTON EXPANSION AND THE QMPI 

Let us initially consider a model in which the term in the action bilinear in the quantum 
field $ is of the form. 



where p = —id and and V are functions of the external fields. The one-loop generating 
functional is hence given by 



In Euclidean space, a representation of F^^^ is provided by Schwinger's proper-time integral 
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r(^) [A„ V] = lr ^Str e-^* (2.3) 
2 Jo t 

{H = \{p — Af + V^), or, upon regulating F*^^) using the ^-function 
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An expansion due to Schwinger 

Stre^+^ = Str [ /^dai 5 (1 - ai) e"^^ + 
1 

/ daidaa 5(l-«i -a2)e°i^5e°2^+ , (2.5) 
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can be used in (p.3|) or ( p.4|) to extract contributions to F*^^) to any order in and l^, 
thereby allowing us to compare one-loop Green's functions. The Str operation, if performed 
in momentum space, results in a loop-momentum integral while integrals over the paremeters 
ai are direct analogues of Feynman parameters integrals. In fact, from (|2.3| ), the expression 
for the A^-point function is identical to that obtained from the Feynman diagram approach 
once the denominators from individual propagators are combined. However, we note that 
the regulated expression for A^-point function obtained by using the Schwinger expansion 
of (|2.5|) in conjunction with the regulated one-loop generating functional of (|2.4| ) cannot 



be obtained by evaluating Feynman diagrams using an action that has been regulated in 
some way; the functional determinant in ( |2.2| ) is what has been regulated and hence the 
designation "operator regularization" is used in 

In thermal field theory, fields are periodic in Euclidean time if they are Bosonic and 
antiperiodic if they are Fermionic; hence time-like momentum variables in these two cases are 
of the form n {2t[ / (3) and (n + 1/2) (2t{ / (3) respectively [IH. A sum over n (the "Matsubara 
sum" ) replaces the usual integral over the time-like component of the loop-momentum. It is 



possible to convert this sum into a contour integral |jTO[ and then to extract the temperature 
dependence (T = 1//5) by considering a forward on-shell elastic scattering amplitude with 
an appropriate statistical weighting factor. To illustrate this, let us consider the particular 



case where in (|2.1|) , = and V = |-F(x), a background scalar field in D dimensions. 



From (2^) then 



, ^tre-K^^+^)*, 
2Jo t 



r(i) [F] 

which to second order in F gives, using 
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Computing the functional trace in (2l) in momentum space as in we find after inte- 
grating over t that 



d^g 



du- 



F{p)F{-p) 
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(2.8) 



In ( p.8| ) we see that upon integrating over u the Feynman integral for the contributions of 
the two-point functions to the one-loop generating functional is recovered, 

1 



(2vr)"g2 ^q^pY 
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In thermal field theory this becomes 
1 



j d^-VE f{p.^-J)f(-p 
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^+{^) {q + P)'+{f{n + n 
The sum over n in ( p.lO| ) can be converted into a contour integral using the formula 
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dz [/ {^z) + f{-^z)] 

f(tz) + fi-tz) 



dz 
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provided / (iz) has no singularities along the imaginary z axis. All the temperature depen- 
dence in ( p.ll|) resides in the second term; hence we find that the temperature dependence 
in ( |2.1CI|) comes from 
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/r\D—l„ 1 flOO+d 
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Performing a partial fraction expansion in ( |2.12| ) we obtain 
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27r)^ ^27riJ-ioo+5 
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dz N(z) 
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(2.13) 
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_2\p + q\\z + po - \p + q\ z + po + \p + q\ ^ 
Since po = 2mm/ [3, we see that (-z) = N {z + po). This permits us to close the contour 
integral in z on the right half plane, and then to use the residue theorem to obtain 



qD-l 



q N{\q\) 
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1^1 

{q + pf {q-pf. 



(2.14) 
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This is the so-called "Barton amplitude" viz an expression the on-shell forward scat- 

tering amplitude with two external fields, weighted by the statistical factor (|g|). This can 
be generalized to handle more external legs, vertices with momentum dependence, and in- 
ternal particles with Fermi statistics 0. In particular, the non-linear four-point interaction 
of photons has been analysed using this approach 

An alternate procedure for examining thermal Green's functions is to use the QMPI. 
In this approach, we use the fact that a path integral can be used to represent the matrix 
element in (jgj) and (g^). 



< x\exp 



-{p-AY + V 



t\y >-- 



' j Dg(r)exp |y dr 



q(0) = y 
per 
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(P is the path ordering) 
where "per" means invariance under the replacement go (t) qo (t) + P. 
As is discussed in ||TO|J^ (p.l5|) can be rewritten as 

1 



q(t)^(xQ + nl3,x) 



l=00 , 

e^'^-P J Vq (r) exp dr 



q' [r)+iA{q{T))-q{T)-V{q{T)) 



(2.16) 



q{0)=y 



with Sn = for Bosons and Sn = mr for Fermions. 

We will restrict ourselves to examining (|2.16| ) to zeroth order in V and second order in 
An (higher orders can be treated in a similar fashion). If 

(2.17) 



(2vr) 2 Jf3An{q{Ti)) =e,exp-{iki-q{Ti)) 



then we have 



q{t) = (xQ+nl3,x) 



Jo Jo (2n) (3 J uo 



(r) --i{T)-q{T) 



q(0) = y 
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where 7 (r) = i 1]^=! [5 (r — r^) ki + 6 {t — ti) e^] provided we keep only terms linear in e^. 

Evaluating this path integral and extracting a high temperature expansion for the two- 
point function from it is outlined in [|rT|]. However, it is not apparent how the results 
one obtains from this approach are at all related to what is extracted from the Schwinger 
expansion of (?^) in conjunction with ( [2.3[ ), and hence the Barton expansion. Indeed, the 
Schwinger expansion of the left side of (|2.15|) to second order in A is 



■^xy ^ X I 



du e-5(i-«)p'* 
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we can also expand ( |2.1(j| ) to second order in A to obtain 

J^y = l^dnJ^'dT^ J2 e*'"/ 'DqiT)expJ^dT[--eir] 



q(0) = y 



(2.20) 



N(ri)-A(g(ri))] [i q {r,) ■ A {q {r,))] 

In order to demonstrate the equivalence of I^y and J^y, we split up the path integral in 
(g) so that 
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J^y = J^dnJ\T2 f; e^^"! d^ld^l Vqir) J (r) jl^g (r) 



q(0)=y i{t2)=^2 <3(ti)=''1 



xe 



(2.21) 



In ( p.21| ) we have replaced the time derivative q {Ti) by (g (xj + Cj 
is vahd in the hmit Cj — > 0. 

The path integrals in ( |2.21| ) can be evaluated using the result 

9(n)=''i 

J Vq (r) exp^ dr (-2^^ (^) - 7 (^) ' 9 (^)) 



g(ri - Ci) /2ei), which 



[27r (n - ra) 
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(2.22) 



1 ri 



drdr'G'(r,r';r2,ri)7(r)-7(r') , 
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where 



G r, r ; r2, n = - r - r - — r + r + + ■ 2.23 

2 2 (ri - r2) n - r2 ri - Ts 



By systematically using ( |2.22| ), (|2.21| ) becomes 
Jxy = dri^ dr2 ^ d6id62 Y.^"' ' ^ 



27r (t - n) 27r (n - r2) 27rr2 



(^1) ^/.2 (^2) 



xexp — - 



{x*-hf , (&1-&2)' ^ (&2-I/)' 



t - Ti Ti - r2 



T2 
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t-n T1-T2 



^1 - ^2 ^ ^2 - y 
n - T2 T2 



(2.24) 
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(we have defined x* = + n(3, x)). Using the Fourier transform 

^ y -ip-x-\p^t 

— e , 



{2Titf'^ J (27r; 
we can rewrite ( p.24|) in the form 

Jo Vo V 2 Mill; 2;y ^27r)^(27r)^(27r)^ 
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t-Ti\ dpi Ti - r2 V dp2 
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Integrating by parts with respect to the derivatives occurring in ( |2.26D leaves us with 
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From the definition of G in ( p.23| ), we see that 



G {t2 + 62, Ti - ei, T2, Ti) = ^\ti - r2 - 61 - 62 | - ^ (ti - r2 - 61 - 62) ^15^ (2.28) 

2 2 Ti — r2 



So that in the hmit ei, 62 ^ 0, ( |2.27| ) reduces to 
|.o-i /"dpi dp2 dp; 



1 ^e^5" e-^^-^*+'^'-M^, (pi -P2) (P2 -P3) 



— / dcTi / d(T2 I n 

4 70 Vo V (27r)'' 



X ^ -S,J-^^^ + (Pi +p.),, (p2 +1>3),,| e-H(--)-^^(--)-^^i-] 



(2.29) 



t 



provided Tj = (Tjt and An is the Fourier transform of An. The Poisson resumation formula 



8 



E /M= E / e^'^™/!/^) (2-30) 

m=—oo m=—oo °° 

can now be used to show that if 6n = nne (e = for Bosons and e = 1 for Fermions) then 
I I |,Lo_M^V (2.31) 



m=— OO m=— OO 



Substitution of (|]3T]) into (|2]2|) leads to 



f^A A /"dpi dp2 dps ^27r / 27r + e/2) 



xi^, (Pi -P2) (P2 -Pa) \-5^.,J-^^^j^ + (Pi +P2)^, (P2 +P3)^J • (2.32) 



It is immediately apparent that (|2.32|) is identical to ( p.l9|) upon insertion of a complete 



set of momentum states into the later, as in PJ^, provided these momentum states have 
Matsubara frequencies associated with their time-like component. One can easily generalize 
this result to show that one can obtain the same result for the matrix element M^y =< 
x\exp — ^{p — A)'^ + V t\y > whether one uses the Schwinger expansion or the quantum 
mechanical path integral for T > 0; consequently the Barton expansion and the expansion 
of the quantum mechanical path integral are equivalent procedures, despite their apparent 
differences. 

Let us now consider how the four-point function can be computed in quantum electro- 
dynamics [|T^] and compare the result with that obtained in 0] where the Barton expansion 
was employed. 

III. THE FOUR POINT FUNCTION 

The one-loop generating function for Green's functions with external photon lines in 
quantum electrodynamics is 
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(3.1) 



If {7m, 7i/} = 2?7;,i. and a 



[7^,7^], then 



.2 ^ 1 



so that in accordance with ( |2.15| ) we must compute 

tr < x\ex'p 



t\v> 



(3.2) 



q{t) = (xQ^n(i,S) 



tr e*"'' F 'Dq (r) exp / dr 

n=-oo <j(a)=!/ ^ ° 



ig^ (r) + (g (r)) ■ g (r) - (g (r)) 



(3.3) 



The term contributing to the path integral in (|3.3| ) that is of fourth order in the background 



field Afj, is 



A = tr y e'"" / Pg(r)e-5/o'i-'?'M / dn H H dr^ H dn 

n=-oo qj)=y ^0 ^0 



4 gi ■ ga ■ ga ■ A3 g4 ■ - - (Fi„^F2,/3 ga ■ A3 g4 ■ A4 + perm.) 
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3^0 

+ Yg [^^'^ai3F2ai3Fz^gFi^^^ + Fi^^ F^^^ F2^g F^^^ + Fi^^F^^^F2^^gF^^^ 

(We have let g^ = g (tj) and A^- = A^ (g(rj)), and "perm." refers to permutations of 1, 2, 
3, 4. The trace identities 

tr = 

tr cT„/3 a^s Faf3 G.,s = -8 tr {F G) 

tr 0-^/3 a^s <yxp Fap G^s Hxp = -32ztr {F G H) 

tr (T,/3 a^s axp a^, F^p G^s Hxp I^u = 64 tr {F G H I -FHIG- F I G H) 
+16[tr(FG) iHI) + {FH) {G I) + (F I) {G H)] 

have also been used.) In order to evaluate the path integral in ( p.4| ), we use the plane wave 
expansion of (|2.17| ). The standard result of ( |2.22| ) and the formula 
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together lead to the contribution of the ^-function to the four- point process to be 

1 
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where = -i {kia(^ip - hpeia), Gij = ||rj - tjI - | (r^ + tj) + nTj/t, G^ 
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Gij. The rescahng Ti = ait allows us to integrate over t; we can then sum over 



n using |13| , P | 

00 
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so that we are dealing with expressions of the form 
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(Here we have e = — 1, 0, 1 and A^ = —ki gij kj, with Qij = jG {cXit, Oji).) Using the integral 
representation [|T5|] 



POO 

Ka iz) = dt cosh (a t) e~ 
Jo 



■z cosh(t) 



(3.9) 
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allows us to evaluate the sum over m in (3^). The Mellin formula 



k+l 



(1 - 2'"^)C (z) ; Cis the contour Re z > 1 - 2e 



(3.10) 



can then be employed; by using it we can integrate over the parameter t used in (|3l9| ) since 
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A final simplification involves using the formula []T3| 
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at last we get 
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With this result we see that (B.6|) becomes 
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(3.14) 



The integrals over at eliminates the first term in ( p.l4| ). Closing the contour C with a semi- 
circle at infinity on the left side of the complex plane allows us to evaluate the integrals over 
z in (p.l4|) by use of the residue theorem. The sum over the contributions of these residues 
gives rise to a power series in = (3, starting at T^. This is consistent with once we 
recall the dependence on (3 in the normalization of the wave function in ( p.l7| ). There is no 
dependence on ln/3 since z(+ (z) {z/2) has only a single pole at 2; = 0, —2, —4, — 6 ■ ■ ■; this 
is consistent with as well. 

In [|rT| it is demonstrated that the vacuum polarization in thermal quantum electro- 
dynamics is transverse; it is anticipated that one can similarly show that the four-point 
function in (|3.14|) is also gauge invariant. 
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IV. DISCUSSION 



In this paper we have demonstrated how thermal effects in field theory can be computed 
perturbatively in the loop expansion by using either a so- called "Barton Expansion" or the 
quantum mechanical path integral. The two approaches are significantly different in their 
technical details; nevertheless we have demonstrated that the results obtained in for the 
four point function at one loop order in thermal QED using the "Barton expansion" are 
consistent with what can be obtained using QMPI. 
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